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Event shapes for Au + Au collisions at 11.4 GeV/c per nucleon were studied over 
nearly the full solid angle with the E877 apparatus. The analysis was performed by 
Fourier expansion of azimuthal distributions of the transverse energy {Ej) measured 
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in different pseudorapidity intervais. For semicentrai coiiisions a pronounced event 
anisotropy is identified beyond that expected due to fluctuations in particle multi- 
plicity. The signal decreases for peripheral and very central collisions. The amplitude 
of the flow signal reaches up to 7% of the mean Ej". 
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Phenomenological analysis [1] of early data and later theoretical predictions [2] 
have led to a picture of nearly complete stopping of baryons in ultra-relativistic nu- 
clear collisions at AGS energies. Recently, these predictions have been experimentally 
confirmed [3-6]. This stopping scenario implies energy and baryon densities of up to 
10 times that of nuclei in their ground state [2,7]. Such large compression factors 
should yield pressure-induced flow effects [8,9], leading to a deflection of nucleons 
away from the high density region. The degree of these flow effects depends on the 
build-up of pressure gradients, on mean field interactions [9] and on the equation of 
state of the hot matter [8, 10]. Their study should provide insight into the proper- 
ties of the dense medium formed in the collision. As a first step towards identifying 
such collective effects we report here the presence of pronounced event anisotropics 
indicating directed sidewards flow of the hot and excited matter formed in Au+Au 
collisions at AGS energy. The results were obtained by analyzing, in different pseudo- 
rapidity (77) intervals, distributions in azimuthal angle (</>) of transverse energy (Et), 
measured with nearly 47r coverage in calorimeters surrounding the target. 

Following predictions based on hydrodynamical considerations [11] flow effects 
in nucleus-nucleus collisions were first observed and systematically studied at beam 
energies of a few hundred MeV/nucleon [12]. The present investigations are focussed 
on much higher energy where a phase transition from hadronic to quark-gluon matter 
is possible. The sensitivity of flow [8, 10] to the equation of state implies that our 
measurements may provide information complementary to the study of single particle 
observables and two-particle correlations generally used in the search for quark matter 
[13]. 

First results for the production of transverse energy in collisions of 11.4 A GeV/c 
197 Au beams with nuclear targets at the BNL AGS accelerator have been reported 
previously [4]. This Au + Au data set, obtained with two calorimeters, the target 
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calorimeter(TCal, -0.5 < rj < 0.8) and the participant calorimeter (PCal, 0.83 < 
i) < 4.5), was used for the flow analysis. Both detectors cover 2tt azimuthally, and 
consist of 832 Nal crystals for TCal, and 512 towers for PCal with individual read- 
out. Events were collected with four levels of PCal Ej triggers. Details of the 
experimental technique and the detector set-up can be found in [4]. The possibility 
to use the calorimetric technique without particle identification for flow measurements 
has been investigated [14] for the system Si+Pb. There it is shown that calorimetric 
measurements sensitive tool to study flow phenomena. 

For the event shape analysis, where the absolute calibration of the transverse 
energy scale is not important, we have not used leakage and response correction as 
done previously [4, 15]. Rather, emphasis was placed on very careful gain equalization 
requiring azimuthal symmetry of the energy deposit averaged over many events. In 
addition, to correct for "dead" channels and other instrumental asymmetries, we 
required that the centroid of the first and second harmonic (for definition see eq. (1) 
below) of the azimuthal Ej distribution be zero on average. Our event shape analysis 
was motivated by a concept for flow analysis tailored to ultrarelativistic collisions [16]. 
A more general approach based on Fourier expansion of the azimuthal distribution[17] 
is employed in our data analysis. 

The E802 collaboration has previously reported [18] a correlation between particle 
production ratios and transverse momenta of projectile fragments (spectators) for Si 
+ Au collisions. The pseudorapidity coverage of our detectors is such that target 
and projectile spectator nucleons account for negligible contributions to the energy 
deposit. Thus the event asymmetries reported here are carried by particles from the 
hot and dense matter in the collision region. 

To study the centrality dependence of a possible flow signal, events were sorted 
into 12 bins, each 20 GeV wide, from 40 to 280 GeV according to the transverse energy 
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deposited in PCal (PCal Ej). Corrections for interactions occurring upstream of the 
target were performed using data taken with an empty target frame. The correction 
is important only at low Et, amounting to 39 % in the 80 - 100 GeV interval, 14% 
in the 120 - 140 GeV interval, and negligible for larger Ej values. After background 
cuts there are about 80,000 events left for the analysis. A possible dependence of the 
event shapes on the angle of the incoming beam was studied using the beam vertex 
detector and found to be negligible. 

For the analysis the pseudorapidity space was subdivided into three intervals: 
backward ( -0.5 < i] b < 0.8 ), middle (0.83 < i] m < 1.85), and forward (1.85 < i] f < 
4.7). Central (pseudo-)rapidity is about 1.7 [4]. The somewhat asymmetric choice of 
intervals is dictated by detector boundaries and statistics considerations. To study the 
event shapes the measured azimuthal distributions were decomposed into normalized 
Fourier components [17]: 



where <f>j is the azimuthal angle of the j detector channel and e 3 T is the transverse 
energy measured in that channel. Fourier components n = 1,2,3,4 were studied. 
For n = 1, v\ describes the radial displacement of the distribution, while ipi is the 



parameter vector and the beam direction. Nonzero values for higher harmonics would 



main point of the following analysis is to determine which of the Fourier components 
v n are significantly different from zero and to study their centrality or Ej dependence. 
Because of the reduced fluctuations in transverse energy production for the Au+Au 
system [4], the transverse energy is closely anti-correlated with the impact parameter. 





angle, relative to the (horizontal) x-axis, of the reaction plane spanned by the impact 



indicate more complex events shapes (n = 2 indicates elliptic components etc.). The 



5 



In Figures l(a,b) we display distributions of the orientation of the reaction plane, 
?/>i, for the backward and forward pseudorapidity intervals for events with PCal Ej 
in the range Ej = 200 - 220 GeV. Both spectra are flat as expected, indicating the 
quality of the symmetrization and gain equalization described above. The sensitivity 
of our data to the reaction plane is evident from Fig. 1(c) where the distribution of the 
angular difference Atjj^ = — tjj{ shows a pronounced forward-backward correlation 
with a peak at AtpY = ir. To quantify the degree of correlation we define the ratio 

R = N(\A^-ir\<ir/2) 

N(\Ail>l f - tt\ > tt/2)' 1 ' 

This ratio R depends strongly on transverse energy, as shown in Figure 2. It reaches 

its maximum for intermediate centrality, but approaches I (implying isotropy) for 

peripheral and the most central collisions. This implies that the best determination 

of the reaction plane is for events of intermediate centrality. Due to larger leakage 

fluctuations for the TCal the best resolution in the reaction plane measurement comes 

from data using only the PCal. Estimates following a procedure described in [17] with 

the shape parameters determined below imply a resolution of 30°. 

In the absence of flow, the central limit theorem ensures a Gaussian distribution 

for d 2 N/dv 2 centered at zero, i.e. 

, v 2 , 



I dN e ^{~7^) l 



<*> (3) 



Nv n dv n a 2 ' ^M\\2(e T ) 2 ' 

Here, N is the number of events, ej is the transverse energy per particle, and M is 
the mean multiplicity of particles in the rj interval for events in a given centrality bin. 
The variance of the Gaussian, <r, is independent of n but is determined by M and 
the €t distribution. The presence of flow will result in a shift by v n of the centroid 
of the d 2 N I dv 2 dt/j n distribution. Assuming that flow does not change the variance a 
one can integrate over ip n to obtain 
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1 dN 

N v v dv, 



n 



= -^exp(- 



2a 2 




n 



) 



(4) 



where I is the modified Bessel function. The distribution specified by eq. (4) will 



shapes is, consequently, the observation of such a minimum in the data. However, a 
minimum at zero is not required for the extraction of v n . 

For each of the three pseudorapidity intervals the measured transverse energy 
data were therefore converted into dN/v n dv n distributions for n = 1, 2, 3 and 4. The 
results are shown for the forward pseudorapidity interval in Figure 3 for five centrality 
bins, ranging from peripheral to very central collisions. The most dramatic sequence 
of shapes is observed for n = 1, where a pronounced minimum at v\ = is observed 
in the intermediate centrality intervals, with maxima at zero for the most central 
and the most peripheral intervals. This demonstrates unambiguously the existence 
of directed, sidewards flow in semi-central Au + Au collisions at the AGS. 

The curves in Fig. 3 are the results of fitting eq. (4) to the data, with 4 free 
parameters to fit, simultaneously, the 4 distributions obtained at each centrality. 
The parameters are v n (n = 1 - 4) and a. The parameter v 3 should vanish for 
symmetric systems. Our results are consistent with this expectation. To reduce the 
number of parameters in the fit we enforced v 3 to be zero. The resulting anisotropy 
parameters v^, v™ and for the three rj intervals are plotted as function of centrality 
in Figure 4 for n = 1, 2, 4. From Figs. 4(a) and (c) one can see that both v\ and v{ 
are nonzero over a significant centrality range, and both reach their maximum values 
at intermediate centrality. The shifts go into opposite direction in the backward 
and forward hemisphere, as demonstrated already in Figs. 1 and 2. These results 
demonstrate that, in the backward and forward region, the azimuthal Ej distribution 
in each event is shifted off center. Since we use normalized Fourier coefficients the 



exhibit a minimum at v. 



n 



if v n > y/2a. A sufficient condition for anisotropic event 
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magnitude of v\ shown in Fig. 4 implies that the amplitude of the directed transverse 
energy flow reaches 7 %. The magnitude of the directed flow gradually decreases 
towards very central or very peripheral collisions. The flow signals do not fully vanish 
even for the most central collisions, as expected because of fluctuations between the 
impact parameter and any global observable. 

We have used the fit parameters uj, v{ and a to calculate the ratio R in eq. (2), 
assuming that there is no other correlation except flow. The results of this calculation 
are plotted as the histogram in Fig. 2 and demonstrate the internal consistency of 
the data. 

In Fig. 4(b), the fit parameters obtained for the rj m interval, i.e. near central 
rapidity, are plotted. As expected in this rapidity interval the n = 1 term corre- 
sponding to directed flow is always close to zero. Small but nonzero v 2 terms appear 
in all rapidity intervals, indicating that the distribution has a slightly elliptical shape 
[16, 17]. Analysis similar to that shown in Fig. 1 shows that the long axis lies in the 
reaction plane. The origin of this elliptic shape is not yet known. It may lie in the 
variation with pseudorapidity of the shift v\ within our finite bin size. Since our mid- 
pseudorapidity bin is not symmetric about central pseudorapidity, we would expect, 
in this case, a nonzero v^, not seen in the data. Elliptic event shapes near mid-rapidity 
may also originate from the pressure gradient during the collision (squeeze-out). In 
that case the long symmetry axis should point perpendicular to the reaction plane, 
not in agreement with the data. Assuming that the elliptic shape is due to the rapid- 
ity variation of the n = 1 term a small perpendicular squeeze-out should be visible as 
a nonzero C 4 . This is actually observed in this pseudorapidity bin although it is very 
small. 

In conclusion, we have demonstrated that using calorimetric methods one can 
determine the orientation of the reaction plane for Au + Au collisions at AGS energy 
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with rather good precision. For the first time we have observed directed flow for 
ultrarelativistic Au + Au collisions. Higher harmonics in the event shape were also 
detected but their origin is unclear, at present. 
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FIGURES 



FIG. 1. Distribution of orientations of the reaction plane for events with (200-220) 
GeV Ej in the PCal using the backward and forward pseudorapidity windows rf (a) and 
r/f (b); (c) reaction plane angle difference distribution. The hatched regions indicate the 
boundaries used to calculate the ratio R (see text). 



FIG. 2. Centrality dependence of the backward-forward correlation R. The his- 
togram was calculated using the anisotropy parameters v\, v{ and a. The quantity 
(Te t = -j^-dEj (top scale) indicates the centrality. For details see text. 



FIG. 3. Distributions of experimental parameters corresponding to harmonics n = 1-4 
(points) and fits using eq (4) for the r/f interval. Each row is labelled by the PCal Ej bin. 
For details see text. 



FIG. 4. Anisotropy parameters v\ (O), % 2 (O) an( i ^4 ( a ) m the three pseudorapidity 
intervals plotted as function of PCal Ej. The top scale is defined in Fig. 2. 
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Fig. 2 
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Fig. 4 



